Manipulating nonequilibrium magnetism through superconductors 
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Electrostatic control of the magnetization of a normal mesoscopic conductor is analyzed in a hybrid 
superconductor-normal-superconductor system. This effect stems from the interplay between the 
non-equilibrium condition in the normal region and the Zeeman splitting of the quasiparticle density 
of states of the superconductor subjected to a static in-plane magnetic field. Unexpected spin- 
dependent effects such as magnetization suppression, diamagnetic-like response of the susceptibility 
as well as spin-polarized current generation are the most remarkable features presented. The impact 
of scattering events is evaluated and let us show that this effect is compatible with realistic material 
properties and fabrication techniques. 

PACS numbers: 73.20.-r, 73.23.-b, 75.75.+a 
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The interplay between out-of-equilibrium transport 
and superconductivity [lj was recently successfully ex- 
ploited in a number of systems in order to implement 
Josephson transistors SSB, 7r junctions 0] and elec- 
tron microrefrigerators [y,Lj], just to mention a few rele- 
vant examples. In this Letter we explore its potential in 
the area of magnetism Q and spintronics and present 
a novel approach to control the magnetization and spin- 
dependent properties of a mesoscopic normal conduc- 
tor. In particular, we show that manipulation of the 
(nonequilibrium) distribution of a normal metal through 
an applied voltage can lead to the control of a number 
of spin-dependent phenomena. The key ingredients are 
superconductor electrodes (with energy gap A) and a 
weak external magnetic field. The resulting Zeeman- 
split superconductor density of states (DOS) was orig- 
inally exploited by Tedrow and Meservey to measure the 
spin-polarization of ferromagnets in the case of Al elec- 
trodes [Tiij . As we shall argue, the interplay between Zee- 
man splitting and nonequilibrium yields dramatic conse- 
quences on quasiparticle dynamics stemming from the 
peculiar shape of the superconductor DOS whose energy 
gap compares well with magnetic fields readily accessible 
experimentally. 

Let us consider the system sketched in Fig. 1. It con- 
sists of two superconducting reservoirs (S) connected by 
a mesoscopic normal metal wire (N) through tunnel con- 
tacts (I) of resistance Ri, The structure is biased at a 
voltage Vc and in the presence of a static in-plane mag- 
netic field H, applied either across the whole structure 
(Fig. 1(a), in the following referred to as a-type setup) 
or localized at the superconductors (Fig. 1(b), b-type 
setup). For the sake of simplicity let us assume a symmet- 
ric structure (a resistance asymmetry would not change 
the overall physical picture). As for the superconductors 
we focus on conventional low critical-temperature thin (< 
10 nm) films. In this case the effect of H on the electron 
spin becomes dominant and, assuming negligible spin- 
orbit interaction the superconductor DOS per spin is 




FIG. 1: Scheme of the structure investigated. An in-plane 
static magnetic field H is applied across the whole SINIS sys- 
tem ((a), a-type setup) or localized at the S electrodes ((b), fa- 
type setup). A finite voltage bias Vc drives the normal metal 
out-of-equilibrium allowing to control its magnetization. The 
N wire is assumed quasi-one-dimensional. 



BCS-like but shifted by the Zeeman energy {Eh 
M*(e) = A/-#|Re[( £ 
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where e is the quasiparticle excitation energy measured 
from the Fermi energy (ep), Afp is the DOS in the nor- 
mal state at Ep (2 spin directions), /is is the Bohr mag- 
neton, and a = ±1 refers to spin parallel(antiparallel) to 
the field. 

At a finite bias Vc, in the presence of H and in the limit 
of negligible inelastic collisions, the steady-state distribu- 
tion functions in the metal wire are spin-dependent and 
are given by [13j 



U(e,V c ,H) = 



(1) 
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where = f {e ± eV c /2), A/J = A^(e + eV c /2), 

Af? = AAf (e - eVc/2), / (e) is the Fermi distribution at 
lattice temperature T and e is the electron charge. Ow- 
ing to the nonequilibrium regime driven by the applied 
electric field, the quasiparticle distributions correspond- 
ing to different spin species behave differently, /+(_) be- 
ing shifted towards lower (higher) energy. The magnetic 
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properties of the N region are entirely determined by 
its (spin-dependent) quasiparticle distribution functions. 
The magnetization density in the wire is indeed given by 



M{Vc,H)=^ B J de{Af?(e)f + (e)-JV N 



(s)f-(e)l (2) 



where AAf(e) = \N N {e F +e + a^i B H) &ndAf N (e) is the 
N region DOS in the absence of magnetic field. The func- 
tion M(Vc,H) vs Vc is displayed in Fig. 2(a,b) for dif- 
ferent magnetic-field values. We assumed a silver (Ag) N 
region (with J\fp = 1.03 x 10 47 J _1 m~ 3 ) at temperature 
T = 0.1 T c , where T c = (1.76 k B )- x A = 1.196 K is the 
critical temperature of bulk aluminum (Al, the material 
forming the S regions) and k B is the Boltzmann constant. 
When H is applied across the whole SINIS structure (a- 
type setup), M. decreases upon increasing Vc starting 
from its equilibrium value M.p a uii = Mb-^f ^ typical 
of a Pauli paramagnet [8j (see Fig. 2(a)). M. shows a 
complete suppression for Vc ^ A/e, i.e. the N region is 
demagnetized. The inset of Fig. 2(a) shows how M{Vc) 
is weakly dependent on the lattice temperature up to 
T = 0.4T C owing to the BCS A(T) dependence together 
with the temperature-induced broadening of /o(e). Con- 
versely, when the magnetic field is localized at the S elec- 
trodes (b-type setup) a negative magnetization is induced 
in the wire (see Fig. 2(b)). Note that M. is antiparallel 
to H . Therefore, the N region behaves as a "diamagnet" 
14]. For eVc > A the wire susceptibility \ (shown in 
Fig. 1(c) at T = 0.1 T c ) reaches the Pauli value but with 
opposite sign x = OM/dH = -n B Mp = ~XPauii- This 
gives rise to a sort of " artificial" Pauli diamagnetism. 

Insight into the physical origin of this 
superconductivity-controlled magnetism can be quali- 
tatively gained by considering the (zero-temperature) 
steady-state DOS diagrams of Fig. 2(a',b'), where the 
normal metal is described by parabolic subbands typical 
of a free electronlike paramagnetic conductor such as 
silver. At equilibrium the occupation of quasiparticle 
states is identical for both spin species leading to 
M. = M pauli and M = for a-type and b-type setups, 
respectively. When Vp ^ 0, electron distributions for 
the two spin populations are characterized by distinct 
chemical potentials [f . Since [i + < fi~~ the occupation 
of spin states antiparallel to the magnetic field is favored 
with respect to the parallel one, owing to the opposite 
energy shift of the superconductor spin-dependent DOS 
in the external magnetic field. This leads to a reduction 
of M. for the a-type setup and to negative magnetization 
for the b-type setup. In particular, at Vc ~ A/e the 
chemical potential separation is <5/.t = /i + — /i~ ~ —2[i B H 
for both setups. This shows a full electrostatic control 
of the mag netization, a unique feature of the present 
system |15j . The superconducting reservoirs are essential 
elements for these effects to be present and replacing 
them by ordinary normal contacts would only always 
lead to a negligible paramagnetic correction to the actual 
M. (second order in the small quantities [i B H/eF and 
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FIG. 2: (color) (a) Magnetization density M vs bias voltage 
Vc at T — 0.1T C for different magnetic fields (H) for a-type 
setup (see Fig. 1(a)). Inset: M vs Vc for different temper- 
atures at H = 0.2 A/ fiB- (b) The same as in (a) for b-type 
setup, (a') Schematic diagrams of the N region density of 
states and quasiparticle occupation both at equilibrium (left) 
and nonequilibrium (right) for a-type setup, (b') The same 
as in (a') for b-type setup, (c) Contour plot of the normalized 
magnetic susceptibility x/XPauli vs Vc and H at T = 0.1 T c 
for b-type setup. 



eV c /e F ). 

The experimental accessibility of this operational prin- 
ciple must be carefully assessed. Electrons in metals ex- 
perience both elastic and inelastic collisions. The lat- 
ter drive the system to equilibrium and can be expected 
to hinder the observation of the phenomena discussed 
here. Our analysis will show a remarkable robustness 
of these effects. At low temperatures (typically below 1 
K) electron-electron scattering and scattering with 
magnetic impurities 111 are the dominant sources of 
inelastic collisions [l&ll9ll2(]| ]. Since Ri is in general large 
compared to wire resistance (Rn = L/ATp e 2 DA, where 
L is the length, A the cross-section and D the wire dif- 
fusion constant), we can assume that f a changes only at 
the interfaces, being essentially constant elsewhere [21j . 
The effect of electron-electron scattering due to direct 
Coulomb interaction on the spin-dependent distributions 
can be accounted for by solving a pair of coupled station- 
ary kinetic equations: 



D 



d 2 h 



dx 2 



n d 2 f-(e) T f , 



(3) 



together with the Kuprianov-Lukichcv boundary condi- 
tions at the NIS interfaces j^. In J3J 2JLi( £ ) is the ne t 
collision rate at energy e, functional of the distributions 
functions f a , defined by 
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FIG. 3: (color) (a) Spin-dependent distribution functions 
/o-(e) vs energy e calculated for three tC co ii values at H = 
0.4A//i B , Vc = 1.8 A/e and T = 0.1 T c . Solid(dashed) 
lines correspond to antiparallel(parallel) spin species, (b) 
fa(s) vs e calculated for various c m values at H — 0.4 A//4b, 
Vc = 1.8 A/e and T = 0.1 T c for a-type setup, (c) The same 
as in (b) for b-type setup, (d) Magnetization density M vs 
Vc at H = 0.4 A/ /i S and T = 0.1 T c for different magnetic 
impurity concentration. Open circles refer to b-type setup, 
filled triangles to a-type setup for c m = 0.001. The latter 
were shifted by — Mp a uli- Data in (b)-(d) were obtained as- 
suming D = 0.02 m 2 s"\ T K = 40 mK and S = \. 



where 



£ n on(e) = [1 - Ue)\ J dJ^U(e w) 

JdE{f+(E + u)[l-f + (E)} + f-(E + oj) [1 - /_ (E)]} 

(5) 



and 



| dE{f + (E)[l -f+(E + u)} + f_(E)[l -/_(£? + w)]} . 

(6) 

In (J5J and @), fc(w) = K ee uj~ 3 ^ 2 according to the the- 
ory of screened Coulomb interaction |23j for a quasi- 
one dimensional wire, where K ee = (nV2Dh 3 ^ 2 J\fp A)^ 1 
E3- By making © dimensionless |2lJ, the 
strength of the electron-electron interaction can then 
be expressed as K, co u = (Ri/RN)(L 2 K ee /D)^/A = 
{L/s/2){Ri/R K )y/K[hD, where R K = h/2e 2 . 

We analyzed quantitatively a realistic Ag/Al SINIS mi- 
crostructure HQ with L = 1 /mi, A = 0.2 x 0.02 /xm 2 , 
and i?/ = 10 3 i7. Figure 3(a) illustrates the effect of 
electron-electron scattering. We solved (jSJ) with H = 
OAA/hb, V c = 1.8 A/e and T = 0.1 T c for several JC coU 



values from negligible (JC co u — 0, blue lines), to mod- 
erate (K, C oii = 1, red lines) and extreme (/C co « = 100, 
black lines) |23l l24| . As expected, electron-electron in- 
teractions have virtually no impact. By increasing the 
strength of Coulomb interaction the quasiparticle distri- 
bution of each spin species relaxes toward spin-dependent 
Fermi functions still characterized by different chemical 
potentials (a similar effect is expected in the presence of 
interaction with the lattice phonons [2(| ) . As a result the 
nonequilibrium magnetization in the normal wire here 
presented is virtually unaffected. 

The situation drastically changes if we assume the pres- 
ence of magnetic impurities in the N region, due to the 
resulting spin-flip processes. Above the Kondo temper- 
ature (Tk), the distribution functions can be calculated 
including, in the right-hand-side of Eqs. (0), an addi- 
tional term derived by generalizing the theory developed 
by Goppert and Grabert in Ref. [23 • It is noteworthy to 
mention that its strength is proportional, apart from the 
electron and magnetic impurity spin coupling constant, 
to the total number of magnetic impurities present within 
the wire (i.e., to the product c m LA, with c m the impurity 
concentration) and The resulting distribution func- 
tions relative to the a-type setup are shown in Fig. 3(b) 
at H = 0.4A//x B , V c = 1.8 A/e and T = 0.1 T c w 120 
mK for various c m values expressed in parts per million 
(ppm). We assumed D = 0.02 m 2 s _1 (typical of high- 
purity Ag), magnetic impurities with spin S = \, and 
Tk = 40 mK (as appropriate, for example, for Mn impu- 
rities in Ag) [2jJ,|30j. By increasing c m , examination of 
the figure immediately shows that spin-dependent distri- 
butions are marginally affected even for impurity concen- 
trations as large as 20 ppm. This shows that in the a-type 
setup the nonequilibrium A4 is relatively insensitive to 
large amounts of magnetic impurities. Figure 3(c) shows 
the /o-(e) calculated for various c m values for the b-type 
setup. In such a case, by contrast, the spin-dependent 
distribution functions tend to merge for much lower val- 
ues of c m thus suppressing the induced magnetization. 
In the presence of a magnetic field across the N region 
(a-type setup) impurity spins tend to polarize yielding a 
suppression of spin-flip relaxation processes for the field 
intensities of interest here |l7L ITsl I27I l28j . This does not 
occur in the b-type setup and makes magnetic impuri- 
ties more effective in mixing spins. The full behavior of 
M{V C ) for b-type setup at T = 0.1 T c and H = 0.4 A/fi B 
is displayed in Fig. 3(d) for several c m values (open cir- 
cles). For comparison, M.(Vc) for a-type setup (filled 
triangles) is shown at low impurity concentration. Wc 
wish to underline the robustness of the induced magne- 
tization, M. being suppressed only for rather large con- 
centrations: the latter can in fact be limited to less than 
0.01 ppm in currently available high-purity metals |30l |. 

These results on the robustness of these effects in re- 
alistic structures make it appropriate to investigate their 
potential for device implementation. An immediate area 
of application is spintronics 0- Let us consider a mi- 
crostructure like that shown in Fig. 4(a), where a nor- 
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FIG. 4: (color) (a) Scheme of the spin-polarized current 
source in the presence of a uniform H. Spin-polarized current 
can be extracted by biasing the middle terminal with Vbias- 
(b) Contour plot of the nonequilibrium current polarization 
Pi vs Vc and H for V bias = 0.1 A/e at T = 0.1 T c . 



mal "probe" terminal is tunnel-coupled to the wire |3lj |. 
Upon voltage biasing with Vbias, the presence of spin- 
dependent distributions in the N region yields a finite 
current polarization Pj defined as 



Pi(V c ,H,V bias ) = 



I/-I-I/4 



where 



Ia{Vc,H, Vbias) = PJ deM*(E)N*{e + eV bl as) 
x[f a (e,V c ,H)- f (e + eV blas )} 



(7) 



(8) 



is the spin-dependent current flowing through the addi- 
tional terminal. Furthermore, p= [2e(A/jp) 2 .Rea:t] and 



R ex t 3> Ri is the probe junction resistance. The cal- 
culated nonequilibrium Pj is displayed in Fig. 4(b) for 
T = 0.1 T c and Vu a s = 0.1 A/e as a function of Vc and 
H . We emphasize that for easily attainable values of Vc 
and H a 100% spin-polarized current consisting of the 
antiparallel spin species can be achieved. A quantita- 
tive estimate for realistic parameters shows that sizable 
spin-polarized currents can be available. For example, 
assuming R ext = 10 4 Q at T = 0.1 T c « 120 mK and 
for an external field of 1.5 T the fully spin-polarized cur- 
rent reaches values up to about 10" 8 A. We stress that 
Pi values largely exceed 50% over a wide region in the 
(Vc,H) plane. Note that at equilibrium by placing the 
wire in external magnetic fields of comparable intensity 
only values of Pj of the order of 10~ 6 -j- 10~ 5 would be 
obtained. 

In conclusion, we have presented a scheme to control 
the magnetic properties of a mcsoscopic metal. Mag- 
netism suppression as well as artificial Pauli diamag- 
netism can be accessed in metal-superconductor mi- 
crostructures thus making available a number of char- 
acteristics of much relevance in light of possible ap- 
plications: (1) Generation of 100% spin-polarized cur- 
rents without invoking the use of magnetic materials; (2) 
full- electrostatic control of magnetization over complex 
nanostructured metallic arrays for enhanced performance 
and optimized device geometries; (3) reduced power dis- 
sipation (10~ 14 -^10~ n W depending on the control volt- 
age) owing to the very small driving currents intrinsic to 
SIN junctions; (4) h igh magnetization switching frequen- 
cies up to 10 11 Hz 21]; (5) ease of fabrication that can 
take advantage of the well-established metal-based tun- 
nel junction technology. 
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